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I I Abstract 

Using the FiNLIE solution of the AdS/CFT Y-system, wc compute the 
anomalous dimension of the Konishi operator in planar J\f=A SYM up to 
eight loops, i.e. up to the leading double wrapping order. At this order a 
non-reducible Euler-Zagier sum, Ci,2,8; appears for the first time. We find 
^ that at all orders in perturbation, every spectral-dependent quantity of the 

K*" Y-system is expressed through multiple Hurwitz zeta functions, hence we 

^ provide a Mathematica package to manipulate these functions, including the 

^ particular case of Euler-Zagier sums. Furthermore, we conjecture that only 

^ Euler-Zagier sums can appear in the answer for the anomalous dimension at 

^ any order in perturbation theory. 

Wc also rcsum the leading transcendentality terms of the anomalous di- 
^ mension at all orders, obtaining a simple result in terms of Bessel functions. 

Finally, we demonstrate that exact Bethe equations should be related to 
^ an absence of poles condition that becomes especially nontrivial at double 

^ wrapping. 
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Interactive feature: If you are reading this article as a pdf file using a 
viewer that supports JavaScript (like Adobe Reader), you can click a sharp 
symbol, for instance this one whenever you encounter it in the text. 
When such symbol is clicked, a pop up window with a Mathematica code or 
example relevant to the context will appear. The Mathematica code requires 
a number of packages from [Tj More examples and further explanations 
how to use packages are given in the notebook usage .nb in [1]. 
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1. Introduction 

During tlie last decade, tliere was a remarkable progress in applying the 
AdS/CFT integrability to solve the planar Af=4: supersymmetric Yang- Mills 
theory (SYM), see [2] for a review. In this context, one of the best studied 
directions is the AdS/CFT spectral problem - the computation of anomalous 
dimensions of gauge invariant local operators or, equivalently, of energies of 
dual string configurations. One typically demonstrates how to solve this 
problem by considering some particular operator or a class of operators; the 
most known example are the so called twist J spin 5* states. The shortest 
member of this family whose anomalous dimension 7 is not protected by 
supersymmetry is TtZD^Z, where D is a light-cone covariant derivative and 
Z is a complex scalar field of J\f=4: SYM. This state known as the Konishi 
operator corresponds to J=2 and S=2. The Konishi operator is interesting 
both from the point of view of N'=^ SYM and its string dual: on the gauge 
side it appears among the leading terms of operator product expansions, 
and on the string side it is among excitations with the lowest energy. From 
the point of view of integrability, the anomalous dimension of the Konishi 
operator is specifically hard to compute perturbatively in the sense that 
the so called wrapping corrections [3] start to contribute earlier than for 
the anomalous dimensions of other states. Explicitly, the leading wrapping 
correction appears at four loops in the perturbative expansion in the 't Hooft 
coupling constant = ^^^Ji" ■ Whereas the anomalous dimensions of "long" 
operators which are free of wrapping corrections (e.g. the ones with J — )■ 00) 
can be studied by solving the algebraic equations of the Beisert-Staudacher 
asymptotic Bethe Ansatz [HE], the presence of wrapping corrections requires 
to solve functional equations instead: the Gromov-Kazakov-Vieira Y-system 
16J. 

Over the past few years, there has been quick progress in the understand- 
ing of how to compute anomalous dimensions of "short" operators (the ones 
that receive wrapping corrections). One already has reasonably precise nu- 
merical values for anomalous dimensions of the Konishi state [71 [8] and of 
various twist J spin S operators [9j. At strong coupling, these dimensions 
were found analytically at two [10] and three [11] loops, with the two- loop 
result matching existing computations fTIi [T^ [T3] from the string theory, 
whereas at weak coupling the analytical answer for J = 2 case was found up 
to five loops: the results [TS, HEl [TTl [lH| for = J = 2 case coincide with per- 
turbative quantum field theory computations that reached four [191 EHl 121] 



3 



and five loop [22] orders, while the analytical continuation of the arbitrary S 
result [23l[2ll|25] to S* = —1 agrees with the prediction from the BFKL equa- 
tion [26]. The Konishi state was recently analysed even in greater detail and 
by now the six [2^ and seven loop [28j results for its anomalous dimension are 
available. There are also improvements in analysing other operators than the 
twist J spin S operators [221 EO], in studying twisted [211 ESI [331 EH ES] and 
g-deformed [Ml ET] versions of the AdS/CFT spectral problem, and in com- 
puting angle-dependent cusp anomalous dimensions and related quantities in 
A/'=4 SYM from boundary thermodynamic Bethe Ansatz [38l [391 SO]- 

Though all these advancements look encouraging, a worrisome sign is that 
a number of important results was obtained by approaches that are stiff for 
improvement and generalizations. 

At strong coupling, using the method of [TO] seems technically infeasible 
beyond two loops; for instance, the three-loops result yjj was obtained us- 
ing two-loop findings and an extra knowledge about the Basso's slope [H]. 
This situation is somehow similar to the computation of the cusp anoma- 
lous dimension by taking a certain limit of the generalized scaling function 
[i2l [131 [B]. While this approach becomes a burden beyond two loops, the 
cusp anomalous dimension can be computed by other means [l5l [16] to any 
desired order, and we might hope to find such means for the dimensions of 
short operators as well. 

At weak coupling, the analytic five-loop result in [IB] was obtained by 
a perturbative solution of the thermodynamic Bethe Ansatz equations, and 
there are conceptual obstacles in generalizing the method of [18\ to higher 
loops [U]. Another approach of [T5l [TBJ is to use Luscher formulae. It is an 
efficient way to account for the so called single wrapping effects. However, 
the Luscher formulae can be derived only for the vacuum state, while their 
generalizations to excited states (in particular, to the Konishi state) is a 
conjecture. At double wrapping orders, Luscher formulae are known and 
were used for the vacuum state of the 7-deformed theory [51], however it is 
not clear how to proceed and to generalize them to excited states [IH]. Hence, 
there is currently a theoretical bound for applying Luscher formulae, which 
is seven loops for the case of the Konishi anomalous dimension and which 
was reached in 

The limitations of these methods and results indicate that some funda- 
mental properties of the AdS/CFT integrable system are still not under- 
stood or were not used in these computations. In particular, the above- 
mentioned five-loop result at weak coupling only partially simplifies the orig- 
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inal way to find the spectrum from tlie tliermodynamic Betlie Ansatz equa- 
tions |19], EOl El], wliereas, in later developments, more simpler structures 
were found behind them: these equations were shown to be derivable from 
the Gromov-Kazakov-Vieira Y-system (or equivalently the Hirota T-system 
on the T-hook), if one supplements it with certain constraints E3] on 
discontinuities of the Y-functions. Then, these constraints were reduced in 
[M] to simple group-theoretical conditions on the T-functions. 

The T-functions obey the Hirota equation, which is a generalization of 
character identities where a dependence on a spectral parameter is added. 
The Hirota equation can be solved by the so called Wronskian solution, a gen- 
eralization of the Weyl character formula, which allows one to parameterize 
the T-functions in terms of a finite set of Q-functions. These properties were 
for a long time known for various integrable models [55], and their equiva- 
lents were discovered for the case of the AdS/CFT integrability in |56l |57] . 
see also [581 EH ED] and references therein for a more generic set up. 

By using these "algebraic" findings, the Wronskian parameterization of 
the T-functions and the group-theoretical constraints on them, a finite set 
of nonlinear integral equations (FiNLIE) which allows one to compute the 
anomalous dimensions of several operators was derived in [51] . By contrast, 
the thermodynamic Bethe Ansatz equations form an infinite set of equations. 

In this paper, following our previous work [27] , we apply the FiNLIE, and 
hence the findings discussed above, to show that the weak coupling expansion 
can be carried to an arbitrary order in perturbation theory. We demonstrate 
this by an explicit computation of the Konishi anomalous dimension up to 
eight loops, a benchmark order at which the double wrapping effects become 
important for the first time. 

The goal of this work was however not only to improve the methods for 
weak coupling expansion. Understanding of the AdS/CFT integrability is 
far from being ultimate. In particular, it is conceivable that, instead of the 
mirror thermodynamic Bethe Ansatz resulting in an unwieldy infinite set of 
equations, there should be a better way to derive the Hirota system and 
constraints on it. In our ongoing research [HI] we are finding new interesting 
features behind the AdS/CFT T- and Y-systems. On the one hand, these 
features should yield a more efficient way to compute anomalous dimensions. 
On the other hand, we hope that they will shed additional light on the funda- 
mental aspects of the AdS/CFT integrability which would also be important 
beyond the spectral problem. It appeared that having explicit analytical 
results is necessary to guideline our research, and that not all of the effects 
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manifest themselves at single wrapping orders. Hence we had to solve the 
AdS/CFT Y-system up to a double wrapping order, and we present such 
solution in this paper. 

Our analysis of the perturbative weak coupling behaviour of the FiNLlE 
showed that all its quantities can be expressed through the so called multiple 
Hurwitz zeta functions. In the following section [2] we acquaint the reader 
with these functions and with basic operations on them. Section [3] describes 
the FiNLIE adjusted to weak coupling expansion and gives some simplified 
examples which show how to use it. All further technical details are avail- 
able online [T] in the Mathematica notebook format. In section |4] we give our 
main result: the eight-loop Konishi anomalous dimension, and in section [5] 
we make some crosschecks by performing an expansion in the inverse powers 
of transcendentality of zeta functions. Finally, we summarize our results in 
the conclusions section, where we also discuss how the exact Bethe equations 
are related to the regularity conditions on a solution of the FiNLIE. The an- 
alyticity properties encoded into these exact Bethe equations are an example 
of a property which acquires new qualitative features at double wrapping. 

To set up terminology, let us note that we define the anomalous dimension 
7 of an operator as 7 = A — Aq, where A is the total conformal dimension 
and Aq is the classical dimension. For the Konishi operator Tt ZD^Z on has 
Aq = 4. The energy E of the string configuration dual to the Konishi state 
is defined with respect to the BMN vacuum TiZ^ so that E = 7Konishi + 2. 
The Konishi operator is a member of the so called Konishi supersymmetric 
multiplet. Another member of this multiplet is X]^=iTi'$?, which is an R- 
symmetry singlet formed from real scalar fields $j of J\f=4: SYM[^ As all the 
operators from the Konishi multiplet, this operator has 7 = 7Konishi5 though 
its classical dimension is different: Aq = 2. 

2. Multiple Hurwitz zeta functions 

In section [3] we will show that all the quantities of the FiNLIE, and hence 
of the Y-system, are expressed at any given order of weak coupling expansion 
through multiple Hurwitz zeta functions. The goal of the present section is 
to define these functions and to list their essential properties. We start with 
a specific case of great importance: Euler-Zagier sums. 



^One defines Z as e.g. Z $1 + i$2- 
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2.1. Euler-Zagier sum^ 

The Euler-Zagier sums, also known as multiple zeta values (MZV), are 
defined as follows 

(1) 



0<rai<n2<...<nfc<oo ^ ^ k 

k 

w = ^ Oj is called the weight, or transcendentality, of the sum, k is called 

i=l 

the depth of the sum. The sum is convergent for > 1, + ak-i > 2, . . . , 

k 

Y2ai > k. In the following we will define a regularization which allows us to 

1=1 

define MZV in the marginally divergent case, i.e when strict inequalities > 
above are weakened to >. We also assume that > 1, otherwise the sum 
can be straightforwardly reduced to sums of lower depth. For instance: 



OO 712 — 1 

C-l,a = X] 5Z = 2 ''^""^ ~ ■ 
712=1 ni=l ^ 

Stuffle algebra. MZVs form a ring over Q. Indeed, consider for example 

n,m \n<m n>m n=m/ 

Clearly, if we consider two arbitrary MZVs, we can repeat this logic to split 
the sum and express their product as a linear combination of other MZVs. 
The equalities that we obtain this way are known as stufHe relations. 

Shuffle algebra. At the same time, there is another inequivalent way to ex- 
press a product of MZVs through a linear combination of other MZVs. For 
this we use the Feynman representation of Euler-Zagier sums: 

Cai+l,a2+l,---,afc+l ~ 



oo>tfc>...>ti>0 



^For a more comprehensive discussion of this subject the reader may consult [62 l l63 l l64] 
and references therein. 
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In the stuffle case we had to spht summation whereas here we will split 
domain of integration, e.^ 



ClC2= / rf/it^tH / dflt2+ Ci/i ((t2-tl)+tl) =2Cl,2 + C2,l, 

Jti,t2 Jti<t2 Jti>t2 

(5) 

where dfx is the measure of integration defined in Q. 

Using this logic, we define the so called shuffle relations. This name comes 
from the fact that we shuffle the ordering of the integration variables in all 
possible ways. In the mathematical literature a different integral representa- 
tion is preferred, the one used to define multiple polylogarithms |65] , however 
the net result is the same. 

Both the shuffle and stuffie products are implemented in our Mathematica 
package zetaf unctions .m [Ij. ^ 

Combinations of shuffle and stuffle products allow one to generate non- 
trivial relations between MZVs. A classical example is that the relations 

C1C2 =0,2 + (2,1 + C3 (stuffie product) , 

C1C2 =2 Ci,2 + (2,1 (shuffle product) 

give the Euler relation 

Cl,2 = Cs . (7) 

Diophantine conjecture. It is conjectured [HS] that all the algebraic relations 
between Euler-Zagier sums are generated by shuffle and stuffie relations. One 
immediate consequence of this conjecture is that MZVs of different weight are 
algebraically independent. This conjecture called the diophantine conjecture 
is not proven. However, one can test it by exploiting integer relation algo- 
rithms in experimental mathematics (see |67] and references therein). These 
algorithms allow performing an efficient and systematic search for relations 
of the type ci(i = 0, where the summation is over a set of multi-indices 
/ and c/ are integers. In particular, these algorithms can show that no re- 
lation exists with coefficients c/-s smaller than a certain magnitude. The 
diophantine conjecture found a solid support using this approach. 



^In this example, both (in the l.h.s.) and C2.1 are marginally divergent. However, 
the formal manipulations that we write disregarding the convergency issue give a correct 
result if one uses the regularization prescription defined in Appendix A.l[ In general, 
the shuffle algebra that we derive this way is valid for a product involving at most one 
marginally divergent MZV in this regularization prescription. 



8 



Based on the diophantine conjecture, one can show [HS] that the number 
My, of independent irreducibk^ MZVs of weight w can be found from the 
following generating relation [68]: 1 — — a;^ = 11 (1 ^ x"^)^^ , which gives 

the following values of My^ when w < 13: 



(8) 



w 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


My, 





1 


1 





1 





1 


1 


1 


1 


2 


2 


3 



In the eight-loop computation that we present here, it is enough to find 
all relations up to transcendentality 13. This is easily done by brute force 
based on shuffle and stuffle relations. In zetaf unctions .m, the result is 
saved into the substitution rules subZetaReduce 7^. We obtain that the 
following MZVs are independent: 



C2 



TT 



(y, for W odd, C2,6 , C2,8 , C2,10 , Cl,2,8 Cl,2,10, C 



1,3,9, 



Cl, 1,2,8 • 

(9) 



As expected, this is consistent with the values of My, tabulated in (|8j). We 
will actually find out that the weak coupling perturbative expansion of the 



Konishi anomalous dimension (63) up to eight loops contains only a subset 
of the irreducible MZVs ([9]), namely it only involves Ci,2,8 and single-indexed 
MZVs. 

If one wants to go to higher weights, one can use more efflcient algorithms 
based on remarkable relations derivable from the shuffle and stuffle algebras 
[691 ITOl ITTl 172] . Assuming diophantine conjecture, all the relations among 
MZVs were worked out in |72j up to weight 22 and, with some restrictions 



on the depth, to higher weights. These results are available online 



2.2. Multiple Hurwitz zeta functions (rj -functions) 

In the evaluation of Feynman integrals, one often encounters polyloga- 
rithms. These functions can be defined as a generalization of Euler-Zagier 
sums which preserve the shuffle but not the stuffle algebra, see e.g. [65] . 



''MZV is called irreducible if it cannot be written as a combination of MZVs of lower 
depth or weight. 



^The marginally divergent case is not included in this statement, see Appendix A.l 
for its treatment. 
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Quite remarkably, the perturbative weak coupling solution of the FiNLIE 
is expressed in terms of functions which are in a sense complementary to 
polylogarithms: they are generalization of Euler-Zagier sums which preserve 
the stuffle but not the shuffle algebra. They are defined as follows 

1 



.U 



E 



0<ni<n2<...<rtfc<cxD 



{u + ini)°-^{u + in2Y 



(10) 



We introduced i = y — 1 and started the sum from so that these functions 
have poles at position 0, — i, — 2i, .... This way, u will coincide with the 
spectral parameter of the AdS/CFT integrable model. This sum is defined 
under the same restrictions on as for the case of Euler-Zagier sums. 



The MZVs are then related to the functions (10) evaluated at the point 



u 



Cai, 



k 



ai,a2,...,ak 



The single-indexed functions (10) differ from the Hurwitz zeta function 
Ca{u) = S„>i(w + a)"" only by a slight change of variable, and they are 
related to derivatives of the polygamma function ip by 

-^('^-i)(-in). 



Va[U) 



a-1] 



(12) 



Hence, the multiple-indexed functions (10) can be called multiple Hurwitz 
zeta functions or multiple polygamma functions. These sums can be also 
recast in terms of infinite cyclotomic harmonic sums defined e.g. in [73]. To 



simplify further presentation, we will often refer to (10) as 77-functions. 



In Appendix A.l which defines marginally divergent 77-f unctions, we 
choose to enforce (12) even when a = 1, which means that we define rii{u) = 
iip{—iu). Relation (11) is also enforced for marginally divergent case, hence 

A number of the FiNLIE functions do 



one aenneS (,1 = 7Euler-Maschcroni- 

depend on this regularization prescription, however it is a good check of our 
computation that the physical quantities, like the anomalous dimension, do 
not depend on the regularization. 

In the following we will denote multi-indices by the capital letters I, J, and 
K. One additional useful property of 77-functions, which holds in particular 
for regularized sums is 



[2] 1 [2] 



where rif'^ =rii{u + i n/2) . (13) 
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?7- functions have an infinite ladder of poles at u = —in, n G In the 
FiNLIE, this structure of poles is natural at weak coupling <^ 1: in this 
limit, the Zhukovsky cuts on the intervals [—2g — in, 2g — in] collapse into 
poles. In the functions entering the FiNLIE, one can in principle expect 
additional poles at shifted Bethe roots, i.e. at u = Uj{g) — in. For the 
Konishi state, there are two Bethe roots, ^1(5') = —U2{g), and Mi(0) = 
At weak coupling, the appearance of such poles would indicate the presence 
of objects like rii{u + Uj{0)) or even of more complicated functions of type 



(16). However, we experimentally observe a fine tuning which results in 
the cancellation of such ladders of poles at shifted Bethe roots, even though 
a finite number of poles may survive for some quantities. Therefore, the 
possible class of functions entering the FiNLIE solution (at weak coupling) 
is strongly constrained, in the sense that the poles of 77-functions correspond 
only to the collapse of Zhukovsky cuts. Further discussion of this question is 
given in section [3^ 

2.3. Integrals involving rj-functions 

The integrals we encounter in this work reduce to a sum of the following 
elementary ones: 

du r_,i 1 



vrjzr^vf, (14) 



00 



— 27ri (u — V 



where //-functions are the complex conjugate of r^-functions, in the sense that 
fji{u) = {r]i{u*))* , where * denotes the complex conjugation. 

An important property of this integral is that it always evaluates in terms 
of r^-functions and MZVs. For example, for Im (f ) > 0: 

du 21 1 [21 6 (^1 (^3 A / \ ^ / \ / \ 

^ Z^^3 J^^V2 =-^--^-,+S^CsV2iv)+C2Vsiv)+Vsd-l 

(15) 

and in the limit f — + i this integral evaluates to — ^Cs • # 

To prove this property, we introduce an intermediate function that will 
be called generalized T^-functionj^ 

" / 1 \ 

^Sa:::./;^H = E Hi— • (i^) 

' ' ' ^-^ ^^\u + inj — Vnl 

0<ni<ra2<---<n(;<oo j=l ^ j j / 



^Generalized 7y-functions are defined by riY^' '^{u) = ^ryj^^ ' "^(u*) 



11 



The integral (14) is evaluated iteratively by taking residues: if the multi- 
index / has the form I = b,I, and if < Im(f ) < 1 , one haQ 

"+00 




— Res 



Res ri\ 

u=v 




u-v) 



[2] 



(17) 



On the one hand, the residue in (17) is expressed in terms of //-functions 
at point u = V (because the derivatives of //-functions are themselves rj- 
functions). On the other hand, the integral term gives rise to two different 
cases: if J = 0, then by closing the integration contour upwards this integral 



{u,0,0,...}\[+21_[-2] 



evaluates to — ReSu=o ^^^^ — ^ — , whereas if / has at least one element, 

it can be written as / = c, /, allowing to move the integration contour again 
to get 



+00 



du 
^27ri 



{vflfl,...} 



h2] 



-[-2] 



du 
^27ri 



\0,vfl,0,...} 
'lb,a,J 



-[-2] 



[u + ioy 



Res 

u=0 



{v,0,0,...} 



[+2] 



-[-2] 



u" 



(18) 



If the initial multi-index / has n elements, the process stops after n iterations, 
and one gets the integral (14) in terms of standard r^-functions evaluated 
either at point u = i (where they are equal to MZVs ([Il|)) or at point u = v, 
and of generalized ?7-functions of the form evaluated at point 

u = i. 



^For simplicity, we assume that the integral is convergent. In practice, sometimes one 
needs to combine few elementary integrands to achieve convergence. 
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As explained in Appendix A. 2 the generalized 77-functions of type 



^{^..,o,o,-u,o,o,...} expressed through standard ?7-functions and //-functions 

evaluated at point v. Moreover, the result of integration is analytic above 
or below the contour of integration and hence cannot contain simultaneously 



//-functions and //-functions. For instance, in example (15) the answer is ana- 
lytic in the upper half-plane and hence it should be expressible only through 
//-functions. [Appendix A. 3 gives an additional class of relations which allow 



one to express //-functions in terms of //-functions (or the opposite) and to 
insure this last property we discussed. This class of relations is also necessary 
to show that terms like the //-functions evaluated at Bethe roots cancel in 
the final expression for the anomalous dimension. 

The integration algorithm and the properties of //-functions described above 



and in Appendix A are implemented in our Mathematica package [T] which 
was used to perform all analytical computations needed to solve the FiNLIE 
at weak coupling. Examples of usage are given in the file usage. nb. 

3. Set of equations for weak coupling expansion 

Our strategy for weak coupling expansion was summarized in |27j^ The 
goal is to determine three quantities p,p2,U from which all the Y-system 
can be straightforwardly restored. As shown in figure 1 in |27], the set of 
equations divides into two iterative processes: on the one hand, there is an 
internal "Y-cycle" which has to run every time we want to compute one 
more order in g"^ for the densities p and p2- On the other hand, there is a 
"wrapping cycle", which computes the function U, and which only needs to 
be run once for every four orders in (/^. 

Compared to the weak coupling expansion [27J, which was performed up 
to six loops, we will make two essential updates. First, we add an equa- 



tion (42) for an auxiliary quantity h which becomes important at 7 loops. 



Second, we will compute the wrapping correction to U using (46) and (44), 
to reach the double- wrapping eight-loop order. By contrast, the asymptotic 
expression for U was sufficient at single-wrapping. 

For the sake of completeness and to provide some technical details, we 
will write down all the necessary equations even though some of them were 



^The reader might find it useful to take a look on ^27 before going into details of this 
section 
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written in f27\ . The reader should however consult [27] or section 2 
of [M] for the notations we use, in particular for the definition of integral 
kernels, Zhukovsky variables, mirror and magic sheets. The equations pre- 
sented below can be used for a straightforward weak coupling expansion to 
an arbitrary order. Although the expansion is technically quite evolved, we 
developed mathematical techniques which allow handling it. 

3.1. Restoring the Y- and T -systems from p, p2 and U 

The T-(Y-)systems involve infinitely many T-(Y-)functions, and an im- 
portant point towards a finite set of equations is to rewrite these functions 
in terms of a finite number of Q-functions [57j, which are themselves param- 
eterized by a few densities [M]- In the finite set of equations (FiNLIE) of 
|54j . the T-functions are parameterized as follows by the functions p, p2, U: 

The T-functions of the right band ({T)j : s > a}) are computed in a 
certain gauge which is denoted by the letter T. They are expressed in terms 
of the function p alone: 

To,. =1 , Tm =s + /C.ip , T2,s =flT^ . (19) 

On the other hand, the T-functions of the upper band ({T^ : a > |s|}) 
are computed in a different gauge denoted by the letter and they are 
parameterized by a handful of g- functions including = qi23/(li- 

or _„[+«] J-a] I [+a]-[-a] [+a] [-a] [+a] [-a] 

-^a,l —1\ q2 + (l2 Ql + 13 + (l4 13 5 

^ _ l+a] -[-a] l+a] -[-a] _ [+a] -[-a] _ [+a\ -[-a] _ [+a] -[-a] _ [+a] -[-a] 

^a,0 — Q'l2 1l2 ' ?34 %4 1l4 1l4 ^23 ^23 1l3 ^24 ^24 ^13 ' 

(20) 

These g-functions are expressed through the functions qi2,qi,q2, and U as 
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follows 



524 =gi2 > 

VM2/ 



(21a) 



gi3g24 - gi4«?23 ^2 - ^2 

gi2 qi2 

13 — T , 94 



(21b) 



9l2 9l2 

The ^-gauge is chosen such that 

(?i =1 , gi2 =(« - Ml - «)(« + Ml + a) , (22) 

where Ui is the exact Bethe root and a is a complex constant, which is 
adjusted to enforce the vanishing ^,o(wi ± i/2) = 0. This constant has 
order a oc an is suppressed by one wrapping (asymptotically, we have 
(gi2)as = Q = — ul), hence it can be found iterativeljj^ Last, the function 
q2 is parameterized by 

q2 = -iu + /C*p2 -}C*Wp^ , (24) 

where we define Wa = ^l^"^?! + Ia'^^qI W = lim VFa for -u G Z (this 
function is then analytically continued to any complex value of u by avoiding 
short cuts) and Wp^ = ^(H^'"*"*^' + H^'""!) on the real axis. We will also use 
W = lim Wa for u E Z later 011^° The term involving W is suppressed in 

a-s>0 I J 

wrapping and can be found recursively. 



a 



^More explicitly, ^ifl{ui + i/2) = implies the following equation: 

(=^1,0 ("1 + 5) - 9i,2("i + i) 9i,2("i)) + 2 i uia'^ + a a {-i + a + 2 ui){i + a + 2 ui) 



-4i(iii - ^)ui 

(23) 

where the first term in the numerator is of order and the remaining ones are of order 
lONote that Wp^ = W = q'^s^^kl^"^ + qlt^ki °' for w e Z but Wp^ ^ W iorueZ. 
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The parameterization (24) is chosen in such a way that, on the one hand, 
the finiteness of the support of p2 (i.e. the condition that p2 = on Zq) is 
equivalent to the condition S/'qi = (which is a consequence of Z4-symmetry 
115]), and that, on the other hand, p2 is of the form ^/Ag'^ — u^fiu), where / 



is an analytic function in the vicinity of the real axi^^ The latter property 
is a consequence of the built-in assumption that Q- and T-functions of the 
AdS / CFT integrable system have branch points of the second order. 

The formulae presented above are nothing but an explicit realization of 
the Wronskian solution of the T-hook [5_7J bisected into semi-infinite bands, 
without taking into consideration any gluing condition on these bands so 
far (i.e. that a global gauge exists in which the Hirota equation is satisfied 
everywhere on the T-hook, including diagonals a = ±s). The parameter- 



ization (19) and (24) was inspired by an approach in [TH [75]; it is just a 
suitable way to represent the Q-functions respecting basic analytical prop- 
erties of the Y-system, with the details about particular state hidden into p 
and p2- Hence, so far we were exploiting merely algebraic properties of the 
Hirota equation, with minor analytical input. Now we come to more physical 
constraints which will determine U, p, and p2 for the Konishi state. 

3.2. Y- cycle 

We use two integral equations for Yi i and ^2,2. First, an equation for the 



product Fi, 1^2,2 is obtained from the analyticity of the ratio 
upper half-plane which implies |5l] the relation 



1 -^1.0 ; 



in the 



log 



1 ^1.0 



X 



K * 2Re 



log 



V 



X 



(25a) 



^^In [54], a slightly different parameterization was used in which p2 was of the form 
Y^4g2 — v?f{u) + {Ag^ — v?)f{u). By contrast, the present convention incorporates the 
/-term into /C * Wpv . 

^^The real part in the convolution was denoted as 776 in [51] . 
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Similarly, from the analyticity of the ratio 



in the upper 



12,2.3^1,0 \Tl,2-^r,l. 

half-plane, one deduces that the following integral equation for holds: 



log 



^1,1 ^2,1 
^2,2 Ti^2 



log 



0,0 



- /C * 2 i Im (^log 

As the asymptotic expressions 

^1,1 



^0,0 



-2 log 



r+ 
'1,1 



(25b) 



2,2 J 



(^a,0) 



)[+a]n[""] 



^2,2 ^^^2 



(P) 



^2 _ y2 



Aas(^il) 



2 as 



-4 



Aas(Ml) - 4 



(26) 



with Aas defined in (60), solve (25a) and (25b), we subtract them from the 
exact quantities in (25) and get equations more suitable for weak coupling 
expansion. 

Let us introduce the notation [F]^ = and define some standard 

combinations 



^1,0 
^0,0 



r = lo£ 



Ti.i 



r* = log 



Ti,i 



(27) 



which will be used throughout the text. 



By subtracting the asymptotic quantities (|26|) in (25), one gets 



^0] _ 



-0] 



log (^^] =H + 2r-+ - + 2r^'-'^ - 2r^-^+'^) 

V^2,2 1 1,2 J r 



(28a) 



(28b) 



The equations (25 ) and (28 ) are valid in the upper half-plane. The latter ones 
are used in derivation of equation (46) on (f/)r- Together with their complex 
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conjugate, they also determine {X2^2l^2,2)r which appears in the exact Bethe 
equations (50). 



To find the densities p and p2, one needs in particular to analytically 
continue (28) to the domain where u is real and belongs to [—2g,2g\. In 
that case, the analytical continuation picks a residue which cancels some 
inhomogeneous terms in the equations and one get 
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iog(yi,iy2,2), = -4i(7v^r^ 

^1,1 '^2,1 



+00— iO 



H 



log 



^2,2 2 



dv 2v 
oc-to -27Ti{2gzy-v^ 



492 



{H + 2r+) 



(29a) 



(29b) 



where we use z = u/2 g, and the weak coupling expansion is to be performed 



in the double scaling regime: g,u with z being finite. In (29), H (resp. 
r~^) implicitly denotes H{v) (resp. r(f + i/2)), and starting from now, for all 
functions that appear inside an integral, if the argument of a function is not 
specified then it is the integration variable. 

Equations (25) on Yi 1, Y2,2 correspond to the physical constraints, and in 



particular one can show their equivalence to the corresponding TBA equa- 
tions [M]- On the other hand, there is a purely algebraic constraint which 
tells us that Yi, 1,^2,2 can be expressed through p,p2, and U (where U is 
suppressed in wrapping and in this sense is known). Explicitly 



2,2 



l + Y^^^ ^:i,2^l,2 ^lfl^l,2 



1 + ^2:2^ 
1 + 



T+ T~ T+ T 

'2,2 '2,2 I '1,1 '1,1 
"72,172,3 70,172,1 



^[l+0]^[-l-0] 

'1,1 '1,1 



(30a) 
(30b) 



One can check that these equations are also satisfied by the asymptotic so- 



lution (26), so we also subtract the asymptotic quantities to prepare for a 



weak coupling expansion. Finally, by excluding Fi i and 1^2,2 from (30) and 



^■^We use the parity ^,s(— u) = S/'a,s{u) ,Ta,s{—u) — Ta,s{u) , which is satisfied for the 
Konishi state, and the reahty of T-functions in the chosen gauges. These properties imply 
that H{~u) = H{u) and r(u) = f{u) = r{—u). 
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(29) and performing some algebra one gets 



ri'+°l+ri-^-°]-r+-r- + log(5'i,o)r 
4i5( Vl 



E 



oo— iO 



27r2 (2^^)2 -t;2 r 



(31a) 



+ rt-^-"' + r+ + r- + log (51 o)r + 2 log ^ ^ 

This is a set of two coupled equations which constitute the Y-cycle. They 
are used to find p and p2- To uniquely solve them, an additional constraint 
= should be imposed. An explicit expansion of these equations is 
worked out in the Mathematica notebook WeakCoupling.nb while here 
we will discuss an instructive simplified example, by putting Wa = U = 
a = H = 0. We have put to zero the quantities that are of order g^, i.e. 
the quantities which are responsible for the wrapping corrections. In this 
approximation we have f^, = f, and the equations (31) reduce to 

, (32a) 

_o -2m{2gzy-v--^^ ^^^^^ 
with 

where q2 = -iu + /C*p2, (g2)as = -iu + i||, c = )C*p and (c)^^ = i are 

defined everywhere in the w-plane outside the cut Zq. The constants q are 
given by Ci = — a^^_. , C2 = but we will only use that q oc 0{g^). 

Next, we define perturbations 5q2 = q2 — (^2)as' Sc = c — (c)as, 5p2 = 
P2 — (P2)as! = P ~ (p)as expand the equations to the linear order in 



f [1+0] _^ f.[ 


-1-0] 


_ f [1-0] _ 


f[-l+0] 


f [1+0] _)_ f.[ 


-1-0] 


+ f[i-oi + 


f[-l+0] 






+ 2 log 
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the perturbations and to the leading nontrivial order in g. To perform the 
expansion, one uses firstly that close to the cut one has 



(^2 



[±o] 



Cl 



gx 



[±o] 



+ 0{g') , 



X 



[±0] 



(34) 



with = =z± Wl - , and 



-1 



-27ri z — z 



(35) 



where we view 6p as a function of ^ = u/2g G [—1,1]. One will find that 
only a certain combination of 5c, namely 5ct+°] — 5c^^^^ = 5p, appears at the 
leading order. For Sq"^^^ one can use the representation which is analogous 
to (35), however it is more instructive to keep 5g2'^°^ as it is. 

Secondly, one has the following perturbative expansion if the argument 
of a function is far from the cut: 



(^2)as(M) 

6q2(u) 



iu + i^ + 0{g'), (c) 
u 
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2 i J_i TT u — 2gw 



u 



= 0{g'), 
M[6p2] + ... 



(36) 



where we defined M[f] = j\^f{z). 

By comparing (35) to (36), we see that, far from the cut, these functions 



are suppressed by one order of g compared to the same functions close to the 
cut. In particular, 6c with its argument far from the cut does not appear in 
the leading order expansion. 

The argument of a function is far from the cut in two situations. The 
first one is when one encounters a function with a shift: /^^^l = f{2gz ± i) = 
/(±i) + . . . . The second one is when we integrate over the argument, as in 



the integral in (32b). The contour of integration can be deformed such that 



it never approaches the branch points u = ±2g. Hence one can apply (36) 



and take the integral in (32b) explicitly. 



By implementing the outlined strategy, one gets at the leading order 

(xM)^ - (x[-°l)^ = X,6p = i A2 M[6p2]{x^+'^ - x[-°]) , (37) 



where Ai = ^^(^^i ^2-4ci-i) ^ 

^ C2(l-2ci) ^ 



ci(2ci C2— 4ci — 1) 
(l+4ci)(ciC2-ci-l) ■ 
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Note that the equation (x^'^^^)'^ 6q^^^^ — {x^^^^y dqif^^ = always has a 
solution 5q2 = ^{x + I) which has good analytic and parity properties (it 
decreases at infinity and it is antisymmetric on the magic sheet). Hence, for 



any constant M one can solve (37) by 



iX2M A f l\ 



X X^ \ X 



where A is adjusted such that M[6p2] = M. This one-parametric ambiguity 
is fixed by recalling that one should impose ^,i(ui) = 0. Indeed, since 
(=^,i)as vanishes at Bethe roots by construction, one should impose that the 
correction to obeys: 

5^i,i(wi) = 5qt + + W,\ = -^M[5p2] + . . . = 0, (39) 

where . . . denotes corrections of higher order in g. These corrections define 
the value of M[6p2], and in particular the magnitude of 6p2. In our simpli- 



fied case Wi = and hence M[5p2] = 0, which means, due to (37), that 
6p2 = 6p = 0. This is only expected because the asymptotic solution should 
be exact when we set to zero all the terms that are responsible for wrap- 
ping. Although the solution we've got is trivial, its derivation allowed us to 



demonstrate the essential features of the perturbative expansion of (31) 



When we add the wrapping corrections, they produce non-zero terms in 



the r.h.s. of (32) which are suppressed in powers of g and induce corrections 



to the asymptotic expressions for p and p2- At the leading order one has 



6p =^Vl - (-162 - 432 Cs + 432 Cs + 504 Ct), 
6p2 =g'^y/l - 



72 656 
- 378 - 864C3 + -Cl + 48CI + — Cl + 624(5 + I68C7 
5 35 

+ Im (-144^3770+ + 144r/0+ - 24737^0+) + (^432 - 432(2 



+144C3 + — C2' + 576CI + — Cl - 1536C5 - 672(7 ) 



(40) 

where 7^0+ = r?a(4| + |)- 



Since Z^"-* and Z^^"-* are even functions analytic in the vicinity of 
the real axis, they are represented in the double scaling regime as finite degree 
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polynomials in at any given order of g, c.f. (40). This property eventually 



follows from (37) supplemented with sources from wrapping terms, and it 



can be used to transform (37) to algebraic equations for the coefficients of 
the polynomials. 

In WeakCoupling.nb pQ, we compute the first three orders of 5p and 
5p2i 5p up to the order g^^ and 5p2 up to the order g^^ . These orders are 
necessary for our computation of the energy up to eight loops. One formally 
needs also the fourth order of 5p2, however it enters only as M[(5p2], because 



of expansion similar to (36), and hence it is found from (39). 



3.3. Wrapping cycle 

An equation on the function f/, which constitutes the wrapping cycle, is 



X 



-I 1,1^^2,2^0,0 



(41) 



where on the r.h.s we outlined two factors encountered previously in the text. 



Both of them are analytic in the upper half-plane, where (41) is defined. 



This equation on U contains an auxiliary function h which also appears 
in the Bethe equations and which is found from: 



log h = —2 log X + Z* log 



-^+(^-1,11^2,2 - 1) 



J-""*" = Ap Y\ cosh [n [u — Uj)]) exp 



(42) 
(43) 



R\[-2g,2<;] 



df 
27 



(coth[7r(-u - v)] + sign[f ]) log [^1^1^2,2] 



Equation (|41|) actually allows finding U only up to a normalisation which is 

(44) 



fixed by the following additional condition 

1 — ^1,1^2,2 



UU 



^0,0^0,0 



P2^W - VFpv 

For this reason the normalization constant A^? is inessential. 
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The asymptotic expressions for h and U are given by 



where - 



5-i(M+i/2) 

CTl(«-i/2) 



--X ^(Ti, 



=Au 



X 



Bi-) 



B 



{+) 



1-D2 



(45) 



n '^BEs(w, ttj) is the BES dressing phase [5]. 



Let us sketch how these asymptotic expressions are derived. To this end, 



/,[+o]/,[-o] 



let us note that equation (42) is equivalent to the Riemann-Hilbert problem 

(i-yi,iy2,2) ^ under the boundary condition h ~ u~^~^f'^ at 



M — oo. By considering h as function of x (which is justified because h 
has only one cut on the physical sheet), and using the periodicity of J-", 



one can derive from this Riemann-Hilbert problem the equation 



1 —V^ 

1,1 2, 



For the asymptotic values (26) of p and Y-s, this is precisely 
on the dressing phase analytically continued using 



the crossing equation 
the trick of |77j. This explains how the expression (45) for (/i)as is obtained. 
To simplify (41) for asymptotic quantities and to find (f/)as, one uses the 

71,1 



following remarkable relation: 



XV X 



(+) 



1 -^1 B 



(-) 



Like for the Y-cycle, we subtract the asymptotic solution from (41 ) to get 



log 



+oo-iO 



— oo— iO 



U 



V + u 



H 



+ 



1 



1 



u — V + i u + V + i 



u+i / 1 4g2 

' 5 V {u+ir I H 



dv 
^2i7r 



(46) 



\og{h)r = Z* 
+ 



log 

dv 



YiiY, 



2,2 



P 



l0g(Fl,l>^2,2) 



2. ((coth(7^(2(7^^; - v)) - coth(7rt;)) log(ri,ir2,2)r)p. 



, (47) 

where the expression log(Fi iF2,2)r should be computed in the double scaling 
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regime u = 2gz, using (29a 



). Note that 



l0g(yi,iy2,2 - l)r = (t^^t) ^Og(Y,^,Y2,2)r + • • • , (48) 
\ 111 2,2 / as 

and the subleading terms of this expansion (denoted by . . .) are not necessary 
for the eight-loop computation in this paper. 

Finally, the equation for the normalization of U ( [44) ) should be trans- 
formed into an equation for the normalization of {U)r, which requires to 



continue (46) to real values of u and to re-expand it in the double scaling 
regime u = 2g z. For the leading wrapping correction of U, this gives the 
following equation for the normalization of U: 

log(f/l2l^[-2])^ = log| ^^^ii^^io 1 +r+ + r- 

+ 2 / ^ — ^ • (49) 



oo— i 



27ri V V 1 + v 



The dependence on z should be the same in the r.h.s and the l.h.s., which is 
used as a nontrivial check for our computations. 

Let us note that log(/i)r = constant ■ + O {g^^)- Since an overall nor- 
malization of h is irrelevant, log(/i)r does not contribute to the computation 
of U at eight loops. However, the terms in log(/i)r of order g^^ and g^"^ are 
needed to find corrections to the Bethe equations, as we will see below. 

S.Ji^. Bethe equations 

The Bethe equations can be written as 

=1 where e*(") ^ - \ ^4^4^ 




e^^-^^ =1 , where e^^^^ = - | ^ ) • (50) 

-'^2,2 '1,2 II I 

The asymptotic function (0)as precisely coincides with the logarithm of the 
Beisert-Staudacher asymptotic Bethe equations [H [5]. The function is 
subject to corrections to its asymptotic value; these corrections are computed 
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from the following equation: 

/ 



log 




y- y- 
-"1,1 2,2 



^2,1 '1,1 



V 



■^2,2 r2„ 



^2,1 7i 1 



n,! ^2,1 

1-2,2 



(51) 



/ 



he square root i n the numerator (resp the denominator) is expressed from 



(29a) (resp from (29b)), and h is expressed from (47), while 7i,i and ^2,1 = 

yl\ — -f— 21 ' ' 

Q'2 + 5'2 + ^2 are found from the Wronskian parameterization discussed in 
section 13. 1[ 

What one really needs is the value of 5(j) aX u = ui which we parameterize 
as follows: 



50(Mi) 



(52) 



We explicitly computed mi,. . ., m4 which are used for our computation of the 
energy up to eight loops. Their expressions can be found in WeakCoupling . nb 
[T] . Using this data, one finds the exact position of the Bethe root if one no- 
tices that, since 0(mi) = 0, one has (50(mi) = 0(mi) — (0)as(wi) = — (0)as(Si + 
5u\). The latter equation is solved perturbatively as follows 



Ml = Ml +1 ( --g 



^10 + 2g'' - 8g 



12 



9 



,14 



(mi + g'^m2 + g'^m^ + g^m4) 



+ 0{g 



(53) 



where ui is the solution of the asymptotic Bethe equation (0)as(iii) = 0. Its 
explicit expression is 



Ml 



1 

7! 



1 



+ Ag'-10g' + 8g%7 + 3C3)+ ^?'(-461 - 240(5 - 144(3 



+ 4 ^^^(1 133 + 252 Cs + 378 Cs + 630 (7) 

- 6 ^^2(7945 + 1556 Cs + 48 C| + 1944 Cs + 2772 (7 + 4704 (9 

+ 24 ^1^(21577 + 4572 (3 + 240 Cs Cs + 24 C| + 4784 Cs 

+5706 C7 + 8064 Cg + 13860 Cii) + • • 



(54) 
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3. 5. Summary of the perturbative expansion 

To conclude, the analytic perturbative solution of the FiNLIE relies on a 
perturbative expansion of the Bethe root ui, the two densities p and p2, and 
the function U. 

In practice, one computes the deviations from the asymptotic expres- 
sions: 6p = p - {p)^g, 6p2 = P2 - (P2)as) using equations (|3T|, and log (U)^ = 
logU/ (f/)as; using equations (46) and (44). These equations contain a num- 



ber of auxihary objects. Firstly, H,r,f^ (defined in (27)), and log(^^o) 



Wa are computed from the Wronskian parameterization given in section 3.1[ 
Secondly, h and log (1^1,11^2,2)1. are found respectively from equations ( [47| a nd 
(29a). These auxiliary objects appear in the source terms of (31 ), (46), (44)^^ 



and they turn out to be suppressed by one order of wrapping, hence their 
contribution to the solution of the FiNLIE at given order is obtained simply 
by knowing the solution of the FiNLIE at lower orders. 

The position of the exact Bethe root constraints the solution through the 
conditions <5^^(mi) = and i^,i('Ui) = 0. Let us stress that the "asymptotic 
solutions" (26), (45) are defined using the exact Bethe root ui, as opposed to 



the solution ui of the asymptotic Beisert-Staudacher Bethe equation. This 
exact Bethe root is found from equation (50). The correction ui — Ui to the 
position of the Bethe root is suppressed by wrapping and, again, it can be 
computed iteratively. 

3.6. Analytical structure of functions 

It might be not immediately clear that the equations presented above can 
be solved analytically at each order of the weak coupling expansion. We are 
going to show that this is indeed so and we will precise the class of functions 
which appear during this expansion. 

Let us start with the integration of densities p and p2, with either a 



Cauchy kernel /C, like in (19) and (24), or a Zhukovsky kernel like in (42). As 



we saw, these densities are \/l — z"^ times functions which are, at each order in 
perturbation theory, polynomials in z. Actually, in the double scaling regime 
(when z = u/2g is fixed while n ^ 1), which one has to use for integrands 
if they are integrated on the finite support [—2g,2g], any quantity has an 
expansion which is either a polynomial or a polynomial times square root (or 



The perturbative expansion of (44), after subtracting the asymptotic solution is ex- 



pressible through these auxiliary quantities, c.f. eight-loop approximations (49 1 and (48). 
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the sum of these two cases). Such expressions can be exphcitly convoluted 
against both the Cauchy and the Zhukovsky kernel. The result of integration 
is either a polynomial in z, if the value of integral is computed in the double 
scaling regime, or a rational function in u, if it is to be computed in the 
ordinary regime, i.e with g ^ keeping u fixed. 

Let us now consider the weak coupling expansion of the quantities in the 
ordinary regime. It is easy to check that when expanding asymptotic quan- 
tities in this regime, one always gets expressions which can be represented 
as linear combinations of terms of the type fjj ^^_^-^a Vj ■ the follow- 
ing such linear combination will be called a standard-type expression. We 
will now discuss two non-trivial operations one encounters during solution 
of the FiNLIE at weak coupling and show that both of them keep us in the 
class of standard-type expression^^ Hence, in this way we will iteratively 
demonstrate that exact quantities are explicitly computable and that they 
are always expressible in terms of standard-type expressions. 

The first operation is integration. Except for the integrals over the finite 
support [—2g,2g], which we discussed already, all the integrals in the weak 
coupling version of the FiNLIE have an integration contour from — oo to -|-oo. 



parallel to the real axis. In section 2.3 and in Appendix A we showed how 



to integrate standard-type expression along such contour, and the result is 
always a standard type expression. It is easy to iteratively verify, starting 
from the asymptotic solution, that the integrands are always of the standard 
type, and probably the only troublesome place is the equation on h, and 
more precisely the integration with coth.{7T{2gw — v)). For this, we can note 
that 

coth.{'K{2gw — v) = —{rif\v — 2gw) + ^(n — 2gw)). (55) 

Hence, this coth is expressed as a linear combination of r^-functions at each 
order of expansion in g"^. Since ry-functions obey the stufHe algebra and 



hence form a ring, we see that at a given order in g"^, multiplying (55) by a 
standard-type expression always gives a standard-type expression. 

The second operation is a semi-infinite summation. One such summa- 



tion is needed to compute U: if / = log[//[/t^l is known from (46) then 



\ogU = Er=o/'^'''- The other semi-infinite sums are used to compute 



gi3,gi4,g24 according to (21a). In both of these cases, the sums are of the 



^^We do not discuss less complicated operations like algebraic manipulations, under 
which this class of functions is clearly stable. 
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type ^fcLo /'^'^^ some function /. One can already note that, in these 
sums, the corresponding functions / are analytic in the upper half-plane, 
hence the sums can involve ?7-functions, but they should be free of any f]. 
Now, if f-s contain only poles at positions u = —ik, then the sum reduces 
to a standard-type expression, e.g 



-jr^' = (56) 



(u + ik) 



However, if there is a pole in another position, a more generic class of func- 



tions (16) may appear. In principle, one could have poles at Bethe roots, but 
we observe a remarkable cancellation of these poles for the physical solutions 
of the Y-system. 

To give an example of how the cancellation comes out, let us study the 
computation of gis at the leading order, i.e. for t/^ = — ^ and qi2 = Q = 



- uf, where Ui = -U2 = 



[2k+l] 




U — ^ — Uj U + ^ + Uj^ 



(57) 



We see that each pole at a position u = uj+i k appears in two successive terms 
of the sum which cancel each other. We checked this cancellation mechanism 
at the first five orders of the perturbative expansion (which was needed for 
computing the eight-loop anomalous dimension). Hence we conclude that at 
least at these orders the following equation, which is behind the cancellation 
mechanism, holds: 



(58) 



u- 

at the zeroes of qi2 which are u = ui + a and u = —ui — a. At the first 



four orders, we have qi2 = Q and (58) is nothing but the asymptotic Bethe 



equation. At least at the fifth order, (58) is still true, though it is no longer 
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the asymptotic Bethe equation. It is not the exact Bethe equation either, 
although it is equivalent to it because it does not follow from the constrains 
^loi'^i) = = and it is satisfied only if ui is the exact Bethe root. 

In order to have the same cancellation of poles in the sums expressing gi4 
and g24, one should also require that 

=92" points u e{ui + a, —Ui — d} . (59) 



In view of (21b), this requirement can be considered gularity condition 



on At first four orders it is equivalent to the equation ^,i(mi) = 0, and 



we also verified perturbatively that (59) holds at the fifth order. 



Under the assumption that (58) and (59) hold at any order of the pertur- 
bation theory, all semi-infinite sums for qij result in standard-type expressions 
only. 



By analyzing how H and f appears in the r.h.s. of equation (46) for 
U, we see that the cancellation of poles in qij propagates to the statement 
that logf//f/t^l only has poles at position u = —ik, hence U is given by an 
expression of standard type. 



Let us stress that we verified (58) and (59) at the first five orders, and 



then predicted the general form of the functions under the assumption that 
they hold at any order. These two equations look very natural, since they 
ensure a more regular structure of g-functions, by cancelling out some ladders 



of poles. If they are not satisfied, then arbitrary generalized ?7-functions (16) 
would appear. This class of functions still forms a ring which insures that at 
most generalized r^-functions are present in the answer at any order; however 
the answer would be considerably more complicated. 

This cancellation of poles at shifted Bethe roots does not mean that rj- 
functions evaluated at Bethe roots never appear. Such numerical constants 
appear in various places, e.g. in the constraints £^^\,{ui) = and ^,i(mi) = 
0. For instance, they are already present in the leading displacement of the 
Bethe root. The statement that we discussed above is merely that there 
are no r^-functions that depend simultaneously on the Bethe root and on the 
spectral parameter. 



Since the computation of energy using (60) reduces to the computation 



of integrals of standard-type expressions, we conclude that if the constraints 



(58) and (59) hold, then the perturbative expansion of the energy is always 
given only by MZV-s and by //-functions evaluated at Bethe roots. However, 
77- functions at Bethe roots are expected to cancel [15] . We observe that they 
cancel indeed, at least up to eight loops. 
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4. Anomalous dimension 

The anomalous dimension 7Konishi = Axonishi — 4 is found from Axonishi = 
Aas(Mi) + A^rap, where 

AM =4 - 8ylm f 4) ' ^wrap = / ^:^^M£L . (60) 

As compared to the prediction from the asymptotic Bethe Ansatz, A^s re- 
ceives the following corrections through the displacement of the Bethe roots: 

v3 k=i 

(61) 

The perturbative solution of the asymptotic Bethe equation gives 

Aa.(Mi) =4 + 12g^- 48 / + 336 / - 12 /(235 + 24 (3) 
+ 12 ^^"(2209 + 360 Cs + 240 (5) 

- 12 (22429 + 4608 Cs + 3672 (5 + 2520 (7) 
+ 24 ^^^(119885 + 29064 C3 + 144 C| + 24156 Cs 

+ 19656 C7 + 14112 Cq) 

- 12 ^^^(2654761 + 742680 Ca + 5760 Ca Cs + 6624 C| + 623904 C5+ 

528552 C7 + 447552 Cg + 332640 Cii) • (62) 

Separately Aas(Mi) and A„i.ap depend on ?7-functions evaluated at Bethe 
roots, however this dependence cancels out in their sum and one gets the 
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result which only involves Euler-Zagier sums: # 

AKonishi = 4 + 12 _ 48 / + 336 g"" + 96 g\-2^ + 6 Cs - 15 Cs) 

- 96 ^i°(-158 - 72 Cs + 54 Q + 90 Cs - 315 Ct) 

- 48 ^^^(160 + 5472 Cs - 3240 Cs Cs + 432 C| 

- 2340 Cs - 1575 C? + 10206 Cg) 
+ 48 ^^^(-44480 + 108960 Cs + 8568 Cs Cs - 40320 Cs C? - 8784 C| 

+ 2592 Cl - 4776 Cs - 20700 C| - 26145 Ct - 17406 Cg + 152460 Cii) 
+ 48 c/^^(1133504 + 263736 C2 Cg - 1739520 Cs - 90720 Cs Cs 

- 129780 Cs C? + 78408 Cs Cs + 483840 Cs Cg + 165312 C| 

- 82080 Cl Cs + 41472 C| + 178200 C4 C? - 409968 Cs + 121176 Cs Ce 
+ 463680 Cs C? + 49680 C| + 455598 Ct + 194328 Cg - 555291 Cii 

- 2208492 Cis - 14256 Ci,2,8) 

+ . ' ' (63) 

We will now discuss some checks of the consistency of this result. In the 
next paragraphs we will discuss the order of magnitude of the answer, and 
then show to what extent this answer can be predicted from the existing 
numeric data that gives the anomalous dimension for various values of the 
coupling. Another check will be given in the next section, where the terms 
having the highest transcendentality are derived. 



The coefficients of the expansion (63) evaluate numerically to 



AKonishi ^ 4 + 0.75000 (4^)2 - 0.18750 {^gf + 0.08203 (4^)^ - 0.05030 (4^)^ 
+0.03578 (4c/)^°-0.02728 (45)^2+0.02175(4c/)^^-0.01791 [^gf^ . 

(64) 

It is expected (see e.g. section 3 of [44j) that the radius of convergency 
of weak coupling expansions in the AdS/CFT integrable system is at most 
1/4, because various Zhukovsky branch points in the w-plane collide when 



= ±i/4 The coefficients in (64) are in agreement with this expectation. 
One can further verify our results numerically using the 5-digit precision 
data for the anomalous dimension [78j which is given in figure [l] A Pad 



^^There might be other types of singularities which depend on quantum numbers. For 
instance, the radius of convergency shrinks to zero when the spin quantum number ap- 
proaches the value 5 = — 1. 
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0.98 
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Figure 1: Numeric data provided by N.Gromov [78] for the anomalous di- 
mension of the Konishi operator at various values of the coupling g. # 

approximant of a function usually has a larger radius of convergency than 
its Taylor series. And indeed, we empirically observe that the following Pad 
approximant 

12 g'+tckg'" 
AKonishi = 4 + (65) 

k=l 

converges, when increasing A, at least for g < 0.6. We took A = 10, expressed 
the coefficients Ci, C2, ... ,Cq through other coefficients so as to reproduce 
exactly the power series expansion up to seven loops, and then we fitted the 
the remaining q-s and di-s against the numerical data for g < 0.6. With this 
procedure, we got a prediction 0.01790 x 4^^ for the eight-loop coefficient of 
the series expansion, which is in perfect agreement with our analytical result 
(|63|PP 



^^The validity of this approach can be tested by using it to compute 6- and 7-loop 
quantities first. 

'^^In fact, one computes rather AKonishi — Aas(ui) than Axonishi- The magnitude of this 
difference is about 13% of Axonishi, hence the non-triviahty of our test is in reproducing 
the first two digits of the difference. 
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5. Expansion in inverse transcendentality 

A good illustration of the approach explained in section |3] is to perform 
an expansion in inverse powers of transcendentality. Though this expan- 
sion is technically more complicated, there are number of benefits as well. 
In particular, the less number of iterations is needed to capture interesting 
higher-wrapping effects. There is also no need to distinguish the ordinary 
and double scaling regimes, and for the few leading orders of expansion, only 
quite little should be known about the position of Bethe roots. 

The strategy of this expansion is the following. We introduce a variable 
r, which will be a bookkeeping variable for the inverse transcendentality 
order (the transcendentality deficit). We assign this variable to the various 
quantities according to the following rule 

Cai,...,a, r-S-«»Cai,...,a, , (66) 

x,x,z ---^ x,x,z, 
(u + r)'", a>0 r""(u + r)~". 

After r is assigned, we perform an expansion in powers of r. In most cases 
T can be treated as an ordinary variable. However, there are exceptions: for 
instance, when one multiplies {u + ri)~"'^ and [u + r2)~°^ with ri 7^ r2, one 
should perform a partial fraction decomposition (i.e. rewrite the product 
as a sum of terms with a single pole), giving rise to a sum of terms with 
different orders in r, in which the leading term is of order ^jj^g 
transcendentality of m in a numerator depends on whether it cancels or not 
some poles. For us, the latter issue is relevant only for ^2 = — iw + • • •, and 
in this case one can assign r"*"^ for u. Another example of such exceptions is 
that the integration parallel to the real axis may increase the leading power 
of r by 1 and to produce an expression with mixed powers of r. 

To compute the leading transcendentality coefficients in the anomalous 
dimension to all orders in g, let us assume that the leading transcendental- 
ity g-functions can be derived solely from the large volume solution. This 
assumption is justified by a careful analysis of equations which can be found 
in Tanscendentality .nb [Ij. 

Under this assumption, we get 

qi2=Q + OiT'), q, = -^r-' + Oir'), = -%r' + 0{r') . (67) 

6 qx 
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The leading transcendentality term for qi^ is then found from 



oo / 1 \ [2n+l] 

,.3 = -2sV0X:(-^) ^O(r^). (68) 

In this sum, the highest transcendentahty term is determined as follows: at 
each order of g, we write ^2q+q- as a sum of poles ; so that the 

sum over n simplifies to rj^^ [u — Vk), where the most transcendental term 
comes from the pole with the largest exponent. This is the pole coming from 
^, and its prefactor is the value of q+q- at m = 0, namely 9 + O^g"^). Here 
it was enough to know the leading term in the position of the Bethe root 
Ml = + . . ., because the next terms are suppressed in g and hence in r. 
We see that we need to know quite little about the position of Bethe roots 
in order to compute the leading transcendentality contributions. 



Let us define, in analogy with (10 
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0<ni<n2<...<nfc<oo ^ ' ^ J • • • \ J 

Note that fji has the prefactor r°. For qij one gets then 

gi3 = - 18/ r^Qr)^ , q^^ =QigT Qi]^ , 524 =2 Qr)^. (70) 

When we expand H in orders of r, we encounter the following two combina- 
tions: 

Q - ^2 , 9„+ ^ JUn+\ Q - ^13 924 + ^ qti + 524^13 

Si -g:pWl3?2 + ^^14 92 + g24J' ^2 Q+Q- ' ^ ' 

with Si ~ and S2 ~ r^. The expansion of H in orders of r gives 

[Q J" Q / 



^^Thc marginally divergent quantities are defined by fji = X](fTinr ^ u+i n ) ^^'^ 
then by stufHe relations. 
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where the second term is a purely double wrapping effect. We see that at 
the leading order in r, the energy is given by the single wrapping term. 
The above analysis allows computing the highest transcendentality term 



to all orders iia g. It is coming solely from A„j.ap in (60) and is given by the 
following explicit integral: 

Alead t.an = 2i X (-36^) / ^= [\--^ + \\- (73) 

J —2711 /i _ V \ a; a; ' 

-oo-iO V «2 \ 

To this end, we can use the Laplace representations 

1 1 



1 _ V 2;"' 



1,2 



/ dte~'^'dtJa{2gt), (74) 

^0 



oo 



Vf=a{-irl dj^JJl^^ (75) 







which are simultaneously valid for > Im [u] > —1. After some algebra, one 
gets 

.432,^r^a,f^V (76) 



^lead tran 







e*-l~^V t 



In Tanscendentality .nb [T] we also performed a partial analysis for sublead- 
ing transcendentality terms, in particular for those coming from the double 



wrapping term in (72), and confirmed the coefficient of the Ci,2,8 term in the 



anomalous dimension at eight loops. 



6. Conclusions and discussion 



In this paper we solved the AdS/CFT Y-system for the Konishi state 
analytically at weak coupling up to the order where double wrapping effects 
first appear. This allowed us to compute the Konishi anomalous dimension 



up to eight loops (63). At this order we observed the appearance of a non- 



reducible Euler-Zagier sum. In section 3.6 we give a well justified prediction 



that the answer will be a linear combination of Euler-Zagier sums at any order 
of perturbation theory. In comparison, a superficial analysis of Feynman 
graphs allows the appearance of non-reducible Euler-Zagier sums j68j , but it 
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also allows other types of numbers which may appear in the answer starting 



Our computation is based on the FiNLIE that was proposed in [54j and 
was first adjusted for the weak couphng expansion in [27]. We presented all 
the integral equations which allow one to straightforwardly perform the weak 
coupling expansion to any desired order. By contrast, the approach based 
on Luscher formulae is inapplicable beyond single-wrapping order^^ 

It appears that at any order of the perturbative expansion, the FiNLIE 
quantities are expressed as linear combinations of products of a multiple 
Hurwitz zeta function times the complex conjugate of a multiple Hurwitz 
zeta function, where the coefficients of the linear combinations are rational 
functions. We showed that all the integrals can be evaluated analytically 
in terms of this basis. We developed Mathematica packages that handle the 
technical details of the computation and which were used to compute the 
leading double wrapping order. It would be possible to apply the FiNLIE 
approach for other weak coupling computations, e.g for the angle-dependent 
cusp anomalous dimension (c.f. |10]) or the BMN vacuum of the 7-deformed 
theory. We also believe that the Mathematica tools which we provide here 
will be useful in a broader context. 

The main obstacle for even higher-loop computations is combinatorial 
growth of the basis of multiple Hurwitz zeta functions which leads to an expo- 
nentially increasing demand of computing resources. We estimate, however, 
that it is technically feasible to reach triple-wrapping orders if a motivation 
arises for this feat. 

An interesting observation is that the algebra of multiple Hurwitz zeta func- 
tions appearing in the AdS/CFT integrability respects the stuffle relations, 
while the algebra of polylogarithms, more typical for perturbative quantum 
field theory, respects the shuffle relations. Hence, in a sense the two algebras 
are complementary which might be a hint for a richer Hopf algebra structure 
of the planar 7V=4 SYM. 

Instead of performing an expansion in the coupling constant, it is also possible 



■^°We thank Matthias Staudacher for clarifications about this point. 
2ipor vacuum state the double- wrapping Luscher formulae are known [34]. It is not 
clear, however, how to generalize these formulae to excited states. 
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to perform an expansion in inverse orders of transcendentality. This approach 
captures some effects from higher loops, and in particular we computed the 
most transcendental terms to all loop orders ( 76 ) in section [s] Probably 
the higher transcendentality terms may also be captured by perturbative 
quantum field theory methods, which would provide an interesting venue for 
comparison. 



There is another interesting physical phenomenon discussed in section 3.6 



The exact Bethe equations can be interpreted as a regularity condition on 
the Y-system. This feature was already questioned in [31], and its first con- 
firmation was obtained in [27], but only for asymptotic quantities. Here we 
demonstrated that the cancellation of poles in the g-functions is ensured even 
if the leading wrapping correction is taken into account. Quite interestingly, 
we observed a cancellation of poles at the shifted zeroes of qi2 and not at the 
shifted positions of the Bethe roots ±Mi. qi2 is equal to the Baxter polyno- 
mial Q only asymptotically, but wrapping effects make it deviate from Q, and 
then its zeroes acquire a positive imaginary part. On the one hand, qi2 7^ Q 
could be an artefact of the computation scheme. On the other hand, Bethe 
roots are known to become complex in e.g. Lie- Yang model at finite volume 
|80j . hence one can speculate on another interesting possibility: the zeroes 
qi2 are an alternative way to parameterize the physical state. This is espe- 
cially appealing because if we assume that these zeroes are determined by 
the regularity conditions, then equation (58) should be exactly satisfied. On 
the other hand, this equation can be equivalently rewritten as an equation 

-.[+1] 4+1] "[+2] 

^ffr^ = - % for - zero of q,2 , (77) 

which takes precisely the same form as the Bethe equations appearing in 
the algebraic Bethe Ansatz solution of spin chains. Hence the AdS/CFT 
integrable system has one more common feature with integrable systems 
solvable by an algebraic Bethe Ansatz, in addition to the group-theoretical 
interpretation of T-functions proposed in 
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Appendix A. Further details about r7-functions 



Appendix A.l. Marginally divergent rj-functions 



Like in the case of MZVs, the sum ( 10 ) is convergent under the conditions 

k 

ak > 1, ttfe + flfc-i > 2, . . . , ^ aj > A;. By contrast, in the case of rji, the 



1=1 



sum is logarithmically divergent, and we define the regularized value of rji by 



(12). We then define the regularized value of Ci by (11), and this definition 

gives Cl = 7Euler-Mascheroni- 

All the marginally divergent ?7-functions have the form ^7ai,...,afc,i,i,...,i, with 
ak > 1 or k = 0, and we define them by requiring that stuffle algebra relations 
are satisfied. One can check that this definition is self-consistent. 

Explicitly, we have 



Vai,...,a^,l — VlVai,...,ak 
_1 

'7ai,...,afe,l,l =2 I VlVai,...,ak 



i=l 



^ ^ Vai,...,ai-i,l,ai,...,ai^ ^ ^ Vai,...,ai-i,ai+l,ai+i,. 
=1 1=1 
k 

1 ~ ^ ^^ Vai,...,ai-i,l,ai,...,ak,l 
i=l 

k 

' ^ ^ '7ai,...,ai_i,ai+l,ai+i,...,afc,l Vai,...,ak,2 



i=l 



(A.l) 

In this way, we recursively show that all marginally divergent functions are 
expressed as polynomials in rji with coefficients which are linear in convergent 
r^-functions. In particular, there is an explicit formula 



(k times) 



1 



iVi + du)% ■ 



(A.2) 



Using these relations and (11), one defines marginally divergent MZVs. 



Appendix A.2. Factorization property 



As discussed in section 2.3, there is a simple recursive way to express in- 



tegrals of the form (14). After a few iterations, any such integral is expressed 



in terms of standard r^-functions evaluated either at point u = i (where they 



are equal to MZVs (11)) or at point u = v, and of generalized ?7-functions 
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of the form ^jj^--'^'^''"''^'^'---^ ^ evaluated at point u = i. For instance, for the 



integral (15), this procedures gives 

r" + 00 



du 



-[-"A 



(u — V] 



[2] 

:V2 



{v,0} 



, -[-2] [2] 
+ ^3 V2 



(A.3) 

In the expression obtained by this method, it is quite important that 
the generalized r^-functions that appear have only one non-trivial shift v. 
This actually allows us to express them in terms of standard ?]-functions and 
r^-functions. For example, one obtains 7^ 



{vfl} 



+ 2 i Ci - ri2,2 - 2 r]^, 



(A.4) 



The proof of this factorization property is done by induction over the depth 



■If 



of the 77- functions. First, it is clear that r]^^ 

u- 

Let us now assume that, when the depth is smaller than n, we know how to 
express terms of standard //-functions and f^-functions. 

u=i 

Then, for an arbitrary multi-index I = a,I with n elements, we have two 
different ways to compute the integral 



du Vj I 



2] 



-27ri (u 



(A.5) 



for < Im(f ) < 1: either we close the integration contour downwards, and 
we immediately obtain that this integral is zero. Or we close it upwards, 
using the iterative procedure from above. Except for the very last residue, at 
each step we will obtain expressions involving ?7-functions with generalized 
shifts of the form -qj^--'^'^''"'^'^'---^ ^ where K has less than n elements. On the 



other hand, if we denote I 



a„_i, . . . , ai (where / = ai, 02, 

/-{... 0,0, d}\ [+21 



the residue at the very last step is ReStt=o 



the vanishing of the integral (A.5) allows us to express r]j 



-{...0,0,v} 
-{...0,0,'!;} 



, a„), then 



Hence 



in terms 



of usual 77-functions. For instance, this argument allows to derive the relation 



-{0,1)} 
V2,2 



= -C24n2iCi4'^+r/g + 24^! 



u=v 



(A.6) 
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Finally, one notes that the functions ^^■■'^''^'^'^'^'■■■^ where the shift v is not 
in the last position are expressed through the functions using the 

stuffle algebra. For instance, we have 

{0,^,0} _„{0,^} _ {0,0,,;} _ {0,0,,,} . A 7^ 

Va,b,c ~^a,b Vc Va,c,b Vc,a,b ^ ' K^' ' ) 

where the period (■ ■ ■ ) denotes r^-functions with less indices. 



Hence, we showed how to express any basic integral (14) in terms of i]- 
and 77-functions, MZVs, and rational functions of u. 



It is clear that if Im (f ) > 0, the answer for the integral (14) should 
be analytic in the upper half-plane. We will now describe relations that 
allow cancehng all the 77-functions using the periodicity property, and these 



relations can be used to make the analyticity of integrals of the form (14) 
manifest. 

Appendix A. 3. Periodicity property 

There is a certain class of relations between f]- and ?7-functions. One 
example of such relation is: 

(4^'^ + ' = -4C2 (4""'^ + ^2) + (^^'^ + ^4) • (A.8) 



To prove the relation (A.8), we notice that the function (772^^' + 772)^ 



iJ2kGZ (u+ikf ^"^ ^ periodic function with period i, which has the Laurent 

series expansion (^72"^^' + = ^ ~ ^ ~^ ^Cl + Olu"^) in the vicinity of zero, 
which is its only singularity lying in the strip {u : |Im(M)| < 1/2}. In the 



r.h.s. of (A.8), the coefficient —4(^2 is chosen in such a way that the r.h.s., 
which is also periodic with poles on iZ, has the same Laurent series in the 
vicinity of zero. Hence, since also both the l.h.s. and the r.h.s. tend to zero 
when u ^ ±00, they should be equal. 

Using the same idea, we will actually show that an arbitrary function f]i 
can be expressed in terms of the functions t]j and of the function fji. To this 
end, let us study the periodic function 

If / is the empty multi- index, we use the convention = 1, hence we have 
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Let us denote by «2, "^3, ■ ■ ■ i^d the elements of /, and define ii = a. Then one 

[2n2] / 1 \ [2na] 



can generalize the relation (A. 10) as follows: 

[2ni 



1 



1 



ni<n2<...<nd 
d d 

E n 

k=0 ni<n2<...<ni^<0<nk+i<...<nd p=l 
fc=0 



fc=0 
d 



1 



[2np] 



(A.ll) 
(A. 12) 



where we remind the convention 
elements b and c, then we have 



(A.13) 

1. For instance, if / has two 



(A. 14) 



Using the periodicity of Vi, we can then find a set of constants jk such 

that Vi = J2k'yk (^Vi''^^ + ^i) • these constants are found by requiring that 
the singular and constant parts of the Laurent series of the r.h.s. and the 
l.h.s. do match. Using the expression (A.13) of Vi, this allows writing 

til vS+uik+2,-M VhAk-2,-,h- We can therefore 



Efc7fc (^?^^' +^i) -E 



iteratively express an arbitrary function f// in terms of functions rij and of 
the function f]i. 
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